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Abstract. We introduce a novel dark matter scenario where the visible sector and the dark 
sector share a common asymmetry. The two sectors are connected through an unstable 
mediator with baryon number one, allowing the standard model baryon asymmetry to be 
shared with dark matter via semi-annihilation. The present-day abundance of dark matter 
is then set by thermal freeze-out of this semi-annihilation process, yielding an asymmetric 
version of the WIMP miracle as well as promising signals for indirect detection experiments. 
As a proof of concept, we find a viable region of parameter space consistent with the observed 
Fermi excess of GeV gamma rays from the galactic center. 
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1 Introduction 

The existence of dark matter (DM) has been confirmed through many independent studies, 
and we now know that this non-baryonic matter comprises around a fifth of the energy 
density of the universe [1-3]. If DM is composed of a single new particle, then the Weakly 
Interacting Massive Particle (WIMP) paradigm is particularly attractive, since the present- 
day abundance of DM can be set by thermal freeze-out in an expanding universe (see [4] 
for a review). Given that the standard model (SM) features a vast and subtle structure, 
though, there may be comparably rich dynamics in the dark sector. Perhaps the dark sector 
features some of the same properties as visible SM baryons, such as accidental symmetries, 
asymmetric abundances, and instability of some of its components. 

In this paper, we propose a new DM scenario where asymmetries in the visible and 
dark sectors are closely connected. Due to light unstable states carrying baryon number, 
asymmetry sharing is efficient down to temperatures below the DM mass, such that the ul¬ 
timate (asymmetric) abundance of DM is set by thermal freeze-out of a semi-annihilation 
process. This can be regarded as a hybrid framework between asymmetric and WIMP sce¬ 
narios, whereby the processes responsible for asymmetry sharing in the early universe can 
potentially produce signals today in indirect detection experiments. 

As a prototypical example, consider a stable DM particle A which carries baryon number 
1/2. This DM particle couples to an unstable particle B with baryon number 1 and to a 
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Figure 1 . Two approaches to sharing an asymmetry between the DM and SM sectors: (a) through 
a higher-dimension operator that becomes inefficient as the universe cools, or (b) though an on-shell 
intermediate particle B which remains in quasi-equilibrium until A freezes out. In this paper, we 
assume that B chemically decouples from the SM before A freezes out (indicated by the dashing), in 
which case freeze-out of semi-annihilation AA —> Bcj> controls the final DM abundance. 


light mediator cj) with baryon number 0. The semi-annihilation [5-10] process 

AA^Bcj) ( 1 . 1 ) 

allows the SM baryon asymmetry to be shared in the early universe, and it also gives rise to 
promising indirect detection signals in the galactic halo today. Note that this process does 
not involve A. since this is an asymmetric DM scenario. As we will see, there is a viable 
region of parameter space that not only yields the desired DM and baryon abundances but 
also yields a GeV gamma ray signal consistent with the galactic center (GC) excess seen in 
Fermi [11-22], 

The presence of asymmetry sharing is familiar from asymmetric DM models (see [23, 24] 
for reviews), but there is a key difference highlighted in Fig. I. 1 In scenarios shown in 
Fig. la, some high-dimension operator or heavy state links the visible and dark sectors 
[30] (see also [31-40]). As the temperature of the universe drops, this interaction becomes 
irrelevant, locking in independent asymmetries for baryons and DM. In our scenario shown 
in Fig. lb, the interactions that link the visible and dark sectors involve light mediators in 
quasi-equilibrium, and therefore stay relevant to energies below the DM mass. As long as 
tub < 2m A , Eq. (1.1) is always kinematically allowed, impacting both thermal freeze-out in 
the early universe and indirect detection today. 

There are three relevant temperatures in this scenario. The decoupling temperature Tp> 
is when the mediator B chemically decouples from SM baryons, setting the asymmetries in the 
dark sector. The effective symmetry imbalance temperature Tj is when the abundance of A 
(and B) departs significantly from the abundance of A (and B ). The freeze-out temperature 
Tp is when the semi-annihilation process AA —> Bcj) freezes out, setting the final A abundance. 
We will focus on the hierarchy 

T d > T : > T f , (1.2) 

which yields an asymmetric version of the WIMP miracle. As we will see, achieving Tj > Tp 
requires 

m A <m B < 2m A - (1.3) 

1 We assume that the initial baryon asymmetry is determined by some unspecified high scale process, in 
either the visible sector, the baryon sector, or both. There are also interesting scenarios where the baryon 
asymmetry is generated via the freeze-out process itself [25-29], but we will not consider that here. 
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With this mass range it is natural (though not required) to consider B as a bound state of 
two A particles. For example, the semi-annihilation process in Eq. (1.1) could correspond to 
the formation of dark bound states [41], dark nuclei [42, 43], or dark atoms [44], though the 
key difference compared to previous work is that the resulting B particle is unstable. 2 

Unlike standard asymmetric DM models, Eq. (1.1) leads to indirect detection signals 
without requiring a residual symmetric component. The resulting indirect detection spectrum 
depends on the decays of B (and to a lesser extent, the decays of </>). As benchmarks, we 
consider two possibilities involving SM quarks, 3 

B —> udd, B — > ebb, (1.4) 

and both of these decay modes can yield a good fit to the Fermi GC excess, albeit with 
different preferred values of and ms- As we will see, the desired phenomenology requires 
mA and ms to have comparable masses of 0(10-100 GeV). 4 

As a proof of concept, we highlight a benchmark scenario that satisfies the following 
three criteria, summarized in Fig. 2. 

• Correct DM abundance : We require Ddm h 2 G [0.115,0.124], within the 95% Planck 
confidence limit of the DM abundance [3]. 5 

• Correct baryon asymmetry: We require % = (rib — nfi)/s G [8.1, 9.5] x 10 -11 , within the 
95% confidence limit of the baryon asymmetry of the universe [46]. 6 

• Plausible fit to the GC excess : We follow the strategy of Ref. [22] to find the best fit 
masses for the GC excess gamma ray spectrum. 

Simultaneously satisfying these criteria turns out to be possible, since much of the intuition 
derived from symmetric WIMPs holds also in this asymmetric scenario. We check that this 
benchmark is consistent with CMB heating bounds from Planck [3] , and we evaluate possible 
direct detection and antiproton flux constraints. We also investigate searches for displaced 
jets at the Large Hadron Collider (LHC), which could be a distinctive signature for the 
proposed scenario in optimistic regions of parameter space. 

The rest of this paper is organized as follows. We introduce our prototype model in 
Sec. 2 and study its cosmological history and asymmetry sharing in Sec. 3, leaving details to 
the appendices. In Sec. 4, we show how the AA —> Bcj) semi-annihilation process followed by 
B decay can match the observed GC excess gamma ray spectrum. We discuss constraints 
and potential collider signals in Sec. 5 and conclude in Sec. 6. 

2 To achieve a more predictive framework, one could consider B to be a bound state of A particles mediated 
by (f> exchange. In this case, rriB would be determined by mj, m^, and the A-cfi coupling, and there would be 
tight relationship between the various interaction cross sections. In order to explore the full parameter space, 
we will not assume that B is a bound state resulting from cj> interactions. Indeed, in dark nuclei scenarios 
[42, 43], B is bound because of (residual) confining interactions and <p is a spectator to the binding process. 

3 These kinds of decay modes also appeared in the asymmetric DM scenario of Ref. [45]. In that model, B 
itself was DM and had late-time decays. Here, A is DM and B decays are prompt on cosmological scales. 

4 The astute reader may notice that the above processes do not conserve fermion number. This is easily 
reconciled with an extended model presented in App. A, yielding equivalent phenomenology. Similarly, <j> can 
correspond to multiple unstable states, specifically two in App. A. 

’For simplicity, we assume a Gaussian distribution to extrapolate the 2<r Planck range from the 1<t limit 
(9dm h 2 £ [0.1175,0.1219]) in Ref. [3]. 

°Here, s is the entropy density. This differs from the more familiar notation for rjb, which is defined with 
respect to the photon density. Today s = 7.05n 7 [46, 47]. 
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Figure 2. Example region of parameter space where we obtain the desired DM abundance [3] (red), 
the desired baryon asymmetry [46] (blue), and a good fit to the GC excess from AA —> B(f> semi¬ 
annihilation followed by B — > udd (purple) or B —> ebb (orange) decays. Here, an d ijb are 

plotted with respect to the true cross section (< taa^b^v), while the GC fit regions are shown with 
respect to the effective cross section (<jAA^B<j>v) e ff scaled by the DM abundance (see Eq. (4.5)). The 
hashed region indicates the CMB heating bounds from Eq. (5.1), evaluated with respect to ( av) e s. 
Not shown are the AMS-02 antiproton bounds from Eq. (5.5), which do constrain this parameter 
space but have large uncertainties. The white star indicates the benchmark parameters in Eq. (2.7). 

2 Prototype of Dark Matter with a Shared Asymmetry 

The analysis in this paper is based on a simplified dark sector consisting of three particles: 
a stable DM species A with baryon number 1/2, an unstable state B with baryon number 1, 
and a light mediator (/> which couples to the SM through, e.g., the hypercharge portal [48-52], 
Higgs portal [53, 54], or axion portal [55, 56]. The reason the A particle is exactly stable is 
that there are no states in the SM with baryon number 1/2. An explicit Lagrangian with 
these properties is presented in App. A, where the state A is replaced with a fermion/boson 
system. For simplicity, we take this ^a/^a system be mass degenerate in our discussion, such 
that the DM dynamics can be captured by an effective single species A. It is also possible to 
split if )a from 4>a to achieve a more varied phenomenology. 

There is considerable freedom in choosing the masses of A and B, though we focus on 
scales relevant for describing the GC excess: 

'm,A,rnB — 0( 10-100 GeV). (2.1) 

The mass of the <f> is assumed to be small compared to the other scales in the theory, ~ 
0(10 MeV-1 GeV), consistent with <f> being a light dark photon that mixes with the SM 
photon or a light scalar that mixes with the Higgs sector. We take the couplings of 4> to SM 
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Figure 3. Leading interactions that determine the freeze out of A: (a) semi-annihilation which is 
also relevant for present-day indirect detection signals, and (b,c) annihilations which ensure depletion 
of the symmetric components of A and B. 


states to be large enough for cj) to stay in thermal equilibrium in the early universe, but small 
enough to avoid direct detection bounds on A (see Sec. 5). 

The key interactions in this scenario are shown in Figs. 3 and 4. The processes in Fig. 3 
control the freeze-out abundance of A, and for simplicity, we assume these processes have 
velocity-independent cross sections. The semi-annihilation process in Fig. 3a has various 
crossed and conjugate channels of relevance: 

AA —> B(f>, AA —> Bcf>, AB —> A(f>, AB —> A4>. (2-2) 

In order to have an asymmetric DM model with a suppressed symmetric component, we 
assume that the annihilation processes in Figs. 3b and 3c, 

AA —> (j)(j), BB —> cjxj ), (2-3) 

are highly efficient. When needed for numerical studies, we take the amplitudes \M a a^ 0 o\ = 
\A4gg^ ( p < p\ = 1, though smaller couplings lead to similar results. * * * * * 7 In the extreme asymmetric 
limit with negligible A and B abundances, the only (thermally-averaged) cross section of 
relevance is ((Taa^b ^) - 8 

The processes in Fig. 4 control the sharing of the SM baryon asymmetry with the dark 
sector. We denote the total asymmetry by r/ to t, which will later be divided into the baryon 
asymmetry rjb and DM asymmetry r)A • Assuming B is a fermion with decay modes given by 
Eq. (1.4), the B decays in Fig. 4a are mediated by the operators 

-^Bu c d c d c , -4 Bc c b c b c , (2.4) 

A 2 A 2 

1 For such large annihilation amplitudes, one might wonder whether <j) exchange could result in additional 

bound state formation between A and B particles (see, e.g., Refs. [57-59]). Such bound states could be avoided 

either by reducing the (ji coupling strength (perhaps resulting in a small residual symmetric DM component) 

or by ensuring that the would-be Bohr radius of the bound state is larger than the 4> Compton wavelength (see 

App. A). Alternatively, the dominant annihilation process could be unrelated to the <j> particle, as expected if 

A and B were part of a larger strongly-coupled dark sector with both stable and unstable components. 

8 When cj) is sufficiently light, the (semi-)annihiliation process can be Sommerfeld enhanced and therefore 
velocity dependent. This adds non-trivial temperature dependence to the evolution of the comoving A and B 
abundances and can change the viable region of parameter space. We neglect the Sommerfeld effect in this 
work, though one generically expects it to boost the present day semi-annihiliation rate compared to the rate 
during DM freezeout. 
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Figure 4. Leading interactions that determine the asymmetry sharing between DM and the SM: (a) 
decay of B to SM quarks, and (b) scattering of B with SM quarks. 


where A ~ 0(1O-1O 3 TeV) is some new physics scale. 9 The B decay width is given paramet¬ 
rically by 


r B - 


128(2vr) 


m B 


(^) 4 


~ (2.2 x 10 7 sec) 


-l 


/ m B \ 5 / 300 TeV \ 4 
V60 GeV/ V A ) 


(2.5) 


This short lifetime ensures that B decays prior to the start of big bang nucleosynthesis (BBN), 
T ~ 10 MeV, corresponding to f ~ lCP 2 sec [46, 61]. The same operators also contribute 
to the 2 —> 2 scattering processes involving B shown in Fig. 4b, which is the mechanism 
responsible for the initial sharing of the asymmetry between the SM and the DM in the early 
universe (see Fig. 5 below). 

As we will see in the next section, there are three temperatures of interest: the B chem¬ 
ical decoupling temperature Tp>, the A/A (and B/B) symmetry imbalance temperature Tj , 
and the A freeze-out temperature Tp. While Tp> depends sensitively on the scale A, as long 
as T B > max{T/, Tp} then the details of B decoupling are irrelevant for the evolution of the 
Boltzmann system. Furthermore, if Tj > Tp, then the details of the A annihilation are rel¬ 
atively unimportant. Assuming the hierarchy Tp> > Tj > Tp, the dominant phenomenology 
of this scenario can be determined by four parameters: 


{mA,m B /mA,( a AA^B<t>v},r]tot}- ( 2 . 6 ) 

Below, we refer to the following benchmark parameters, which yields the desired behavior 
specified in the introduction with the B —> ebb decay channel: 

mq = 60 GeV, m B /mA = l, {GAA-^B<i>v) = 3.3 x 10 -26 cm 3 /sec, ??tot = 9-4 x 10 -11 . (2.7) 

Note that this benchmark cross section value is comparable to canonical WIMP scenarios. 
The proximity of mq and m B is because the GC excess fit in Sec. 4 prefers a higher Lorentz 
boost factor for B . 10 As shown in Sec. 3.4, we get reasonable DM phenomenology whenever 
m A < m B < 2m a- 

!J If B is a scalar, then the decay must include an additional fermion, perhaps a SM lepton or a supersym¬ 
metric particle (see, e.g., Ref. [60]). If B is a fermionic bound state of a fermion/boson %!)a/4>a pair, then one 
should regard Eq. (2.4) as effectively a dimension 7 operator of the form il)A<t>AU c d c d c , such that the effective 
cutoff scale A would be lower than our benchmark value. 

ll) In addition, as shown in Fig. 7b, taking the ratio ms/uiA to be in the vicinity of 1 gives fine control over 
the DM abundance for fixed semi-annihilation cross section, making it easier to simultaneously fit the GC 
excess and achieve the desired DM abundance. 
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One might be surprised that vtot i n this benchmark is comparable to the present-day 
baryon asymmetry r] b ~ 9 x 10 —11 , such that DM only takes around 1/20 of the asymmetry. 
This is unlike usual asymmetric DM scenarios where the total asymmetry r/ to t is roughly 
twice the baryon asymmetry r/ b [23, 24]. The main reason for the mismatch is that m A is 
around ten times heavier than the standard 5 GeV benchmark for asymmetric DM (i.e. the 
mass that yields the right DM abundance when t/dm — Vb)- Also note that the concrete 
model in App. A has both a boson and fermion species of A, and A has a baryon number 
1/2 (see Eq. (3.1) below), further affecting the expected mass relation. 


3 Early Universe Cosmology 


3.1 Initial Conditions 

We start from a primordial asymmetry generated through an unspecified mechanism (see 
[23, 24] for a review of options). The total asymmetry rj to t is then conserved through the 
thermal history of the universe by the processes in Figs. 3 and 4. We use the notation 
Y t = m/s, where n* is the number density of species i and s is the entropy density. We can 
express the total asymmetry as 

Vt ot =VA + VB + Vb , r < A = \ ( Ya ~ Y a) > Vb =Y b - Y-g, (3.1) 


where the factor of 1/2 accounts for the baryon number of A and r]b is the baryon asymmetry 
(baryons minus anti-baryons). For ease of discussion, we refer to DM as a single particle A. 
Because B is a fermion, however, our explicit model in App. A requires A to be a degenerate 
fermion-boson system, a fact which is reflected in the equations below. 

At high temperatures, the SM and DM asymmetries are related since the operators 
in Eq. (2.4) (or their ultraviolet completions) ensure that the two sectors are in chemical 
equilibrium. A full analysis of the chemical potentials is presented in App. B. In equilibrium 
at a temperature T, the asymmetries are: 


eq 

Va 


hA ( fi{m A /T) + fb(m A /T)) 


Vtot h A ( f{(m A /T) + fb(m A /T )) + h B ft(m B /T) + h b / f ( 0) ’ 

Vb_ = _ h B f{(m B /T) _ 

Vtot h A ( f{(m A /T) + fb(m A /T)) + h B f { (m B /T) + h b / f (0) ’ 

V? = _ hbffjO) _ 

Vtot h A ( fi(m A /T) + fb(m A /T)) + h B fi{m B /T) + h b / f (0) ’ 


(3.2) 

(3.3) 

(3.4) 


where for the explicit model in App. A, {h A , h B , h b } = {1/4, 1, 45/29 ~ 1.6}. (The factor 
of 45/29 is familiar from ordinary asymmetric DM scenarios, see e.g. [62]). The function 
/ is defined in Eq. (B.10); it has the asymptotic behavior f{(x) -A 1/6 for fermions and 
/b(x) —>• 1/3 for bosons as x —> 0 (early times) and f(x) —> 0 as x —> oo (late times). 


3.2 Chemical Decoupling of B 

The equilibrium conditions in Eqs. (3.2)-(3.4) hold as long as the interactions between the 
SM and the B particles in Fig. 4 are active, allowing efficient sharing of the asymmetry. 
When both the Bq —> qq scattering and B -A qqq decay processes go out of equilibrium, then 
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Figure 5. Comparison of the scattering rate of Bq —» qq and the decay rate of B —> qqq with the 
Hubble expansion rate. By Eq. (3.5), the decoupling temperature To occurs when the scattering rate 
falls below the expansion rate, which in turn sets the initial conditions for the subsequent Boltzmann 
evolution of A particles. Chemical equilibrium between B particles and the SM is restored when the 
B decay rate becomes relevant, though in our benchmark studies, this happens after A particles have 
already frozen out. 


the B particles chemically decouple from the SM. To estimate the decoupling temperature 
To, we compare the rate of B scattering/decay to the Hubble expansion 

max{n e ^(T D )exp(-iJ, q /T D )((Tv)Bq-^qq,TB} - H(T d ), (3.5) 

where the thermally-averaged rate for Bq —» qq scattering is calculated in App. C, (av)Bq-^qq ~ 
T 2 / A 4 in Eq. (C.14), the number density of quarks ~ T 3 is given in Eq. (C.9), the chem¬ 
ical potential of quarks g q is defined in Eq. (B.18), and the B decay width T# is given in 
Eq. (2.5). The Hubble parameter is 

H(T) = 1.66^ (3.6) 

where g * is the number of relativistic degrees of freedom at the temperature T [47]. Because 
B interacts with the SM via the contact operators in Eq. (2.4), the scattering process are 
more relevant at early times (high temperatures) while the decay process is more relevant at 
late times (low temperatures). 11 

In Fig. 5, we compare the B scattering/decay rates with the Hubble expansion rate for 
various choices of A, taking the benchmark dark sector parameters in Eq. (2.7). We plot 
temperatures in terms of x = mq/T. For A < 50 TeV, decoupling never happens since B 
and the SM remain in chemical equilibrium throughout the early history of the universe. 
This tends to erode the DM asymmetry r/q, since for T <C mq, /(mq/T) —> 0 in Eq. (3.2). 
Eventually, the semi-annihilation process A A —>• B<f> freezes out, which stops the depletion 
of rjA■ While this could lead to a potentially interesting phenomenology, in this paper we 

11 In addition to processes involving on-shell B particles, there are processes mediated by off-shell B particles 
that can be relevant for chemical decoupling: AA —> qqq and A + q —> Aqq. While these are the same order 
in A as the processes in Fig. 4, they have a subdominant effect due to their 2 —> 3 phase space suppression. 
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(a) (b) 

Figure 6. Evolution of the comoving abundances of the {A,B} system, using (a) the benchmark 
parameters from Eq. (2.7) and (b) a more bound-state-like spectrum with ms/ixiA = 1.9. The 
dashed curves indicate the equilibrium distributions, and the dotted curves indicate the anti-particle 
distributions. The arrow shows the approximate location of the semi-annihilation freezeout yF- In 
(b), it is coincidental that the decay of B starts around the same time as yF- 


focus on larger values of A where B decouples prior to A freeze out. Taking the benchmark 
in Eq. (2.7) as an example, A = 300 TeV yields To ~ 55 GeV, which corresponds to x ~ 1.1 
in Fig. 5. 

When the dark sector decouples from the SM at To, the asymmetries can be estimated 
by evaluating Eqs. (3.2)-(3.4) at To- If To is much higher than rriA,B-, then fi{mA,B/To) -A- 
1/6 and /b(m^/Tp) —>• 1/3, and the asymmetries at the time of decoupling are: 


Vb,D 
m ot 


Vb,D 


0.47, 


= 0.54, 


VA,D 

Vtot 


VA,d 


0.23, 


= 0.17, 


mot mot 

which is not so different from the To —> oo limit. 


VB,D 

-A 0.30. 

Vtot 


actual 

values are 


= 0.29, 

Vtot 



(3.7) 


(3.8) 


3.3 Thermal Freeze-out of A 

Once B chemically decouples from the SM, the dark sector follows standard Boltzmann 
evolution. The full Boltzmann system is described in App. D, which includes the impact of 
having initial asymmetries for A and B. There are two relevant dimensionless time variables, 


x = 


m A 


V = 


mA - 2 m B 


(3.9) 


where the numerator of y was chosen to match the approximate kinetic energy available in 
the AA -A B(f> process. Note that x > y whenever mo < 2 mA- 














In Fig. 6 a, we show the numerical evolution of the various Yx for the benchmark scenario 
in Eq. (2.7). As x increases, the abundances of A, A, B , and B track the equilibrium 
distributions until the A A —> (jxj) and BB -A 0</> annihilation processes freeze out at x ~ 20 (= 
xa) and x — 20 m^/mg (= xb), respectively . 12 At that point, the A and B abundances 
rapidly decrease due to the asymmetry. Assuming mg> tua as in Eq. (1.3), then the effective 
symmetry imbalance temperature is set by xa , yielding Tj ~ tua/ 20. Switching to the y 
variable, the AA —» Bcf) semi-annihilation process freezes out around yp — 20, and Tp > Tp 
since xa > IJf- Eventually, the B particles decay, transferring their asymmetry back to SM 
baryons. We are left with a relic DM abundance of A and a baryon asymmetry rj b , which for 
this benchmark, match the observed values in our universe. 

We now can understand why it is important to take Tp) > Tp. If B were in chemical 
equilibrium with baryons all the way down to Tp then the asymmetry in A would be very 
suppressed. Taking Eq. (3.2) with xa — 20, we find tja — 10 -18 , which is far too small to 
obtain a reasonable DM abundance. For this reason, we need to have B decoupling happen 
before a large A/A asymmetry develops. That said, even though B is chemically decoupled 
from SM baryons at Tp), one can still think of the dark sector as being in quasi-equilibrium 
with baryons out to Tp, since any B particles produced will eventually decay to SM baryons. 
In this way, the asymmetry in the dark sector effectively leaks into the baryon sector via the 
A A -A Bcf) process, and this leakage stops only when semi-annihilation freezes out. 

Perhaps less obvious is why it is important to take Tj > Tp, which requires mg > tua 
since Tp is set by the lower freezeout temperature between AA —> (fxj) and BB -A <p<p. The 
reason is that if there were a large abundance of B at later times, then the AB —> A(j) process 
would still be active during A freezeout, severely depleting the DM abundance . 13 We will 
show this effect numerically in Fig. 7b below. 

In the extreme asymmetric limit and assuming Tp Tj > Tp, we can gain an analytic 
understanding of the DM and baryon abundances. In this limit, the final A abundance is 
negligible, so | Ya{ oo) = 77 , 4 ( 00 ) at late times. After all remaining B particles have decayed, 
the observed baryon asymmetry today is 


?7b(oo) = 77 tot - \ y a { 00). 


Similarly, the present-day DM relic abundance is 

rriASoh 2 


Ddm h = 


PC 


Ya(oo), 


(3.10) 


(3.11) 


where so = 2891 cm ~ 3 is the entropy density today, and pc = 5.15 x 10 ~ 6 GeV/cm 3 is the 
critical density. Thus, in this extreme asymmetric limit, finding the baryon asymmetry and 
the DM abundance reduces to finding 1 ) 4 ( 00 ). 

Moreover, in this limit, the value of Ya( 00 ) can be effectively determined by considering 
just the AA -A Bcf) process. As shown in App. D, the relevant Boltzmann equation is: 

dy A A, , ( v2 ( \ ( y a 1 ) 2 \ , 010 4 

= -^{CTAA^B'PV) ^ - (j7 tot - 7 lb - -Y A J -yeq~J > ( 3 ' 12 ) 

12 For a more accurate determination of the freezeout temperature, one should refer to the full asymmetric 
Xf equation given in Refs. [63, 64]. The large annihilation cross sections considered in this paper have an 
0(1) impact on xa and xb- 

13 When A is composed of multiple states, for example a fermion ipA and a scalar <f>A as in App. A, the 
requirement ms > iua becomes rriB > max{m^,m^ A }. This ensures that the processes ipAB -a 4>\<l> and 
4>aB -A iPa4> are frozen out earlier than the semi-annihilation process 'ipA'pA —> B(j>. 
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where A = s/H, the Hubble scale H(T) is evaluated at T = via — and the equilibrium 

densities are defined in Eq. (D.l). We used Eq. (3.1) to replace Yg with rj to t — Vb — \Ya- 
Similar to the standard WIMP case [47], albeit with the replacement x -A y, Ya stops 
following the equilibrium distribution after freeze-out, and the Boltzmann suppressed term 
approaches zero. Using the freeze-out approximation, the solution to Eq. (3.12) is 


Ya{oo) 


1 A(cr 

Ya(vf) Vf ) 


-i 


(3.13) 


Note that the (Y^ q ) 2 /Y^ q term in Eq. (3.12) scales like e 2y , compared to the WIMP case 
with (Yf) 2 ~ e~ 2x , which explains why freezeout happens at y F — 20 instead of xf — 20. 

We can gain further insights into 1 / 4 ( 00 ) by considering limiting cases. In the limit of 
a large AA —> B<f> cross section, the 1/Ya{uf) factor drops out, and we recover the familiar 
WIMP approximation 

2 mASoh 2 m A s 0 h 2 y F f . 

S2dm h = -Pa(oo) « - — -r• (3.14) 

PC Pc A(<t aa^b<i>v) 


Just as for standard WIMPs, y F has a logarithmic dependence on cross section [47]. In the 
limit of a small AA —> Bcf) cross section, the A and B particles are effectively decoupled 
after freeze-out, so the A abundance is set simply by Ya(uf ), or more accurately, the initial 
asymmetry at decoupling A Ya( oo) — t)a,d • For our benchmark scenario, we end up some¬ 
where in between these two extremes, which is a hybrid of standard WIMP-like behavior (i.e. 
AA —> B<f> freezeout) and asymmetric DM behavior (i.e. initial asymmetries). 


3.4 Results for Dark Matter Abundance 

We now present numerical results for the DM abundance and compare to the analytic ap¬ 
proximation in Eq. (3.13). In Fig. 7a, we show how the DM abundance changes with the 
AA —> Bcj) cross section. For large values of ((jaa^-b^v) , Ddm asymptotes to the standard 
thermal freeze-out expectation (albeit with x —> y compared to ordinary WIMPs). For very 
small values of (a aa^b^v), the A abundance saturates at the initial asymmetry. 

In Fig. 7b, we show how the DM abundance changes as a function of 7715 / 771 , 4 ■ For 
rriA < ms < 2mA, ^dm decreases slowly when decreasing the ms/mA ratio. Once mg< m^, 
there is a dramatic drop in which arises because the AB -A Acj) process is still active and 
relevant for determining 1 ( 4 ( 00 ). It is important to note that when mg > mA, the equilibrium 
term (Y([ q ) 2 /Y^ q decreases sharply, faster than the other terms in the Boltzmann system, and 
therefore setting it to zero after y F is a valid approximation. This does not hold for mg < itia, 
and that is the reason for the difference in behavior between the two regimes. We are mainly 
interested in the extreme asymmetric limit where A and B decouple from the evolution of 
Ya (i.e. Tj > Tp), which is why we focus on the parameter space mA < mg < 2771 , 4 . Again, 
for this mass range it is tempting to interpret B as a bound state of two As, though we will 
not restrict ourselves to that interpretation in this paper. 

In Fig. 8, we show slices through the { 771 ^, 7715 / 771 ^ 4 , {(TAA-^B(j>v), i/tot} parameter space. 
The star marks the benchmark parameters in Eq. (2.7). The red regions correspond to the 
observed DM abundance (Ddm^ 2 E [0.115,0.124]) and the blue regions correspond to the 
observed baryon asymmetry (r/g £ [8.1, 9.5] x 10~ n ). In Fig. 8a, the desired DM abundance 
is achieved with a AA -A Bcj) cross section familiar from the standard WIMP case, with only 
a small variation with rjtot- I n Fig. 8b, the required value of mA also has a weak dependence 
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(a) (b) 

Figure 7. (a) DM abundance as a function of the semi-annihilation cross section (< jaa^B 4 >v)■ The 
curved dashed line corresponds to the approximation in Eq. (3.13), while the two straight dashed lines 
correspond to the extreme freeze-out and asymmetric limits, respectively, (b) DM abundance as a 
function of m b /'tn a ■ Note the dramatic fall off of Ddm h 2 when tub < tua■ In both plots, unchanged 
parameters are fixed to the benchmark values from Eq. (2.7), and the arrow indicates the benchmark 
value for the a; axis. 

on 77 tot, due to the mild logarithmic dependence of yp on via, just as in the WIMP case. Note 
that the resulting DM abundance is somewhat correlated with the total asymmetry, since 
the thermal cross section benchmark is at the transition between the asymmetric regime and 
the WIMP regime, as shown in Fig. 7a. 

Fixing r/totj w e can highlight the mass and cross section dependence. In Fig. 8c, we 
see the presence of two different regimes in agreement with Eq. (3.13) and Fig. 7a. In the 
low cross section limit, the DM abundance is dominated by the initial asymmetry, which is 
indicated by the vertical behavior at rriA ~ 15 GeV. At higher cross sections, rriA has to 
rise linearly with the cross section in order to obtain the correct abundance using Eq. (3.14). 
In Fig. 8d, we also see two regimes which follow from Fig. 7b. For mg/m^ > 1, the DM 
abundance band has only a mild dependence on vtia- As uib/tha — > 1, the number density 
drops dramatically for fixed cross section, so the value of via has to increase to compensate. 

4 Fitting the Galactic Center Excess 

A number of studies have used the Fermi public data to identify a potential excess of gamma 
rays coming from the GC region [11-22], The origin of this excess is as-of-yet unknown, but 
a tantalizing possibility is that this a signature of DM annihilation, though astrophysical 
explanations are also plausible [16, 18, 19, 65-68]. The DM interpretation is bolstered by the 
fact that the needed annihilation rate is consistent with that of a thermal WIMP, though there 
is recent evidence that the GC excess is better fit by a population of unresolved point sources 
[69, 70]. In any case, typical asymmetric DM models do not predict this kind of indirect 
detection signal, though, unless there is a residual symmetric component (see e.g. [71, 72]). 
Here, however, the semi-annihilation process AA —> Bcj) followed by B decaying to hadrons 
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Figure 8. Regions which satisfy the DM abundance and baryon asymmetry requirements: b 2 £ 

[0.115,0.124] [3] (red) and r)b £ [8.1, 9.5] x 10~ n [46] (blue). The star symbol corresponds to the 
benchmark in Eq. (2.7), and from that benchmark value, we sweep (a) (aAA^B 4 >v) versus % 0 t, (b) 
mA versus ??totj (c) mA versus (ctaa^B 4 >v}, and (d) mA versus tob/toa• The dotted (solid) lines 
indicate the lower (upper) experimental limits. The sweep in (c) corresponds to the plot shown in the 
introduction, Fig. 2. 


can give rise to an interesting gamma ray signal. In this way, the GC excess could be 
connected to the process of asymmetry sharing in the early universe. 

The decays of B are prompt on cosmological time scales, so one can think of the AA —> 
Bcj) process as being a one-step cascade decay [51, 52, 55, 73] (see also [74-77]) where the 
B subsequently decays to three quarks. Recall from Eq. (1.4) that our benchmark decay 
modes are B —> udd and B — > ebb, though other flavor combinations are equally plausible. 
These quarks hadronize, and the resulting gamma ray spectrum comes primarily from neutral 
pions which decay via 7 r° —> 77 . The decays of cj) are model dependent and may contribute 
to the gamma ray signal as well. For concreteness, we assume that cf) dominantly decays to 
^ + /r - , 14 which only results in a small contribution to the gamma ray spectrum from final 

14 This is consistent with (j> being a Higgs portal scalar with — 250 MeV. Of course, for larger <j> masses, 
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state radiation (FSR). For simplicity, our study ignores these possible FSR photons as well 
as photons from inverse Compton scattering or bremsstrahlung. 

To determine the gamma ray spectrum in the B rest frame, we use Pythia 8.185 [78] to 
construct a color singlet three quark final state, uniformly filling out the allowed three-body 
phase space. 15 We then use the default hadronization and decay model in Pythia to obtain 
the gamma ray spectrum dN/dE Test after letting all unstable hadrons decay. To go to the 
AA —> B(j) rest frame, we boost the B particle by the gamma factor 


{2mA) 2 ~ + rn 2 B 

4mAmB 



(4.1) 


taking = 0 for simplicity. The resulting gamma ray spectrum is given by (see, e.g., [43]) 

m = rZ/W-m d£ r est dN 

d E 2 / 3 'y JE/{' y(l+/3)) R'rest dE Ies t 

The produced flux of gamma rays as seen on earth is 

d $7 _ T sun d N T ! \ ( a o\ 

dfid E ~ 8tt m\ d E Jnorm (aAA ^ v )^- ( 4 - 3 ) 

The J factor is the integral along the line of sight of the DM density. We adopt the normal¬ 
ization of Ref. [22] where the region of interest (ROI) is |(| < 20° and 2° < |6| < 20°: 


J„„™ = = 2.06 x 10 23 GeV 2 /cm 5 . (4.4) 

J ROI 1 

The effective cross section accounts for the possibility that A comprises only a fraction of the 
total DM density: 


(aAA-+B<t>v)es 


AA-^B^v) 


A A V 


(4.5) 


Written this way, the resulting gamma ray spectrum depends on three parameters 


{m A , m B /m A , (<JAA-^B<j>v) eS }, (4.6) 

and the choice of B decay channel. 

To fit the GC excess, we use the procedure outlined in Ref. [22], Adopting the same 
notation as Ref. [77], the chi-squared for a given parameter point is 

-E^ 

\ O.E j, model (lE j, data 

where CT 1 is the inverse covariance matrix, obtained from Ref. [22]. We show example fits 
(reasonably close to the best ones) of the photon spectrum in Fig. 9, for both the B -A- udd 

the (/> -A 7r°7r° channel opens up, which yields an additional source of prompt gamma rays. We take the Higgs 
portal with mn^ ~ 250 MeV as a benchmark when discussing CMB bounds in Eq. (5.3) and direct detection 
bounds in Eq. (5.6), though strictly speaking, for such low <j> masses, one should also account for additional 
^-mediated bound states of A and B. In the extended model of App. A, we use axion-portal-like couplings, 
where tj> is replaced by a scalar/pseudoscalar pair, in which case there is more flexibility to raise the </> mass. 

15 A more accurate analysis would take into account the matrix element of the B -A qqq decays, but we use 
flat three-body phase space to remain agnostic about the Lorentz structure of the decay operator. 



x 2 = 


V 


,d N 


dN 


E l _ | (J~ 

\ d E i, model 4 B n a 4c / b 


d E i, data 


- 14 - 













Figure 9. Example fits to the GC excess from AA —» Bcj) semi-annihilation with (a) B —> udd 
decay and (b) B —>• ebb decay, following the analysis of Ref. [22]. The black error bars are statistical 
uncertainties while the red band represents correlated systematic uncertainties. The parameters of the 
model are tua, tub/uia , and (aAA^B 4 >v) e s, which makes the number of degrees of freedom 24—3 = 21. 
The effective cross section (av) e g is defined by Eq. (4.5). The solid (dashed) lines are gamma ray 
spectra after (before) boosting from the B rest frame to the AA rest frame. 


and B —> ebb channels. Because the photons from heavy quark decays are softer than for 
light quarks, obtaining a good fit for the B —> ebb channel requires a higher mass than for 
the B —> udd channel. This is consistent with the observation for standard WIMP scenarios 
that b quark final states require higher DM masses than t lepton final states [20]. Because 
we are considering parameter points for which m,A — uig, the effect of the boost is mild, 
pushing the peak of the (energy-squared-normalized) photon spectrum to slightly higher 
values. Consistently with previous work, the best fitting (effective) cross section is close to 
the expected WIMP thermal cross section. 

In Fig. 10, we show the parameter regions that give the best fit to the GC excess in the 
B —> udd (purple) and B —> ebb (orange) channels. In each plane, we first find the best fit y 2 , 
and then show 1, 2, and 3 standard deviation contours for A\' 2 . In Fig. 10a, we show the tua 
versus (cr aa^B 4 > v ) eS plane, leaving fixed, showing that the best fit regions tend to 

have WIMP-like cross sections. In Fig. 10b, we show the jnq versus ms/fnA plane, leaving 
(aAA^B^eS fixed. The value of ms /tua determines the boost of B, and can be used to fine 
tune the gamma ray spectrum. Larger boosts (smaller mass ratios) are somewhat preferred 
by the fit. 

In Fig. 10a, we have superimposed the CMB limits discussed below in Eq. (5.1). These 
limits are largely independent of mg and do not constrain the parameters in Fig. 10b. As 
noted in Refs. [3, 79, 80], models that tend to fit the GC excess are in some tension with the 
CMB limits, though there is still viable parameter space. The CMB limits are proportional 
to a parameter f e g- that quantifies the fraction of energy that goes into electrons and photons. 
We compute this parameter below in Eq. (5.3) and find = 0.17 if we only account for 
the B decay products, though we use the more conservative value f e g = 0.24 which assumes 
4 > -» ■ 
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Figure 10. Fits of the GC excess for the B —> udd (purple) and B —> ebb (orange) channels. 
The best fit regions correspond to regions of 1, 2, and 3 sigma off the best fit y 2 . The star symbol 
corresponds to the benchmark Eq. (2.7). In (a), we fix the mass ratio tub/bia and then vary the 
DM mass via and effective cross section (aAA^B<pv)efi, while in (b) we fix the effective cross section 
(to different values in the two channels) and vary the DM mass and mass ratio. We do not show the 
CMB limits in (b) as those only constrain lower masses. Note that the GC best fit point does face 
potential AMS-02 antiproton bounds (see Eq. (5.5)). The sweep in (a) corresponds to Fig. 2. 


In Fig. 11, we show the same parameter space as Fig. 10, but now letting the third 
parameter from Eq. (4.6) float to give the best fit. In Fig. 11a, we note that floating the 
ratio ms/ruA extends the best fit to a wider range of values for both via and (av ) e g but still 
within the vicinity of the benchmark in Eq. (2.7). In Fig. lib, we see that the best fit for 
each operator is determined by a band in niA, and the change in the mass ratio and can be 
compensated somewhat by a change in the cross section. 

Finally, returning to Fig. 2 from the introduction, we combine the analysis of the DM 
and baryon abundances in Fig. 8c and the GC best fit regions in Fig. 11a. Note that the 
abundances are given with respect to the actual cross section {(Jaa-^b^) , while the GC fits 
and CMB bounds are with respect to the effective cross section (aAA^B<f>v}eff- We made this 
hybrid choice to avoiding display a pathological region of phase space where one obtains a 
good fit to the GC excess with a small cross section but overabundant DM. As advertised, 
the benchmark parameters in Eq. (2.7) yield a consistent cosmology and a plausible fit to 
the GC excess. 

5 Additional Constraints and Signals 

Having seen that we can achieve a viable asymmetric DM scenario with intriguing indirect 
detection signals, we briefly discuss possible additional constraints and signals. 

• CMB heating bounds. The process AA -A Bcj) can occur in the early universe, even 
after thermal freeze-out. This residual semi-annihilation is constrained by limits on 


16 - 











(a v) e ff (10 26 cm 3 /sec) 



Figure 11. Same as Fig. 10, but floating the mass ratio mg/mj in (a) and the effective cross section 
(a A A^B<i>v) e tt in (b). 


the power injected into the CMB through ionization [81]. The power in this case is 
parameterized by: 

feG^ AA^Bd>v)eG 

PCMB = ---, (5.1) 

m A 

where f e g is the efficiency factor and the effective cross section is defined in Eq. (4.5). We 
can divide f e s into contributions from the B decay products and the (model-dependent) 
<p decay products. 

/eff = /Jr + /iff. (5.2) 

Following the analysis of Ref. [82], the efficiency factor from species X is: 


(5.3) 


where we read off the values of /* ff , i e {e^q} from Ref. [82] and estimate /f ff ~ 
0.2(/fg +e + /J ff ) following Ref. [83]. For the B —» ebb decay, we find /Jj = 0.17. For 
the </> decay, we assume that the dominant decay mode is /i + /r _ , leading to = 0.07. 

Current constraints from Planck [3] are shown in Figs. 10a and 11a above (for Figs. 10b 
and lib, the bounds are outside of the plotted region as they constrain rn A < 2 GeV). 
Note that the power injected depends directly on {gaa-aB^)^, and therefore directly 
impacts the GC excess best fit regions. For a fixed effective cross section, the CMB 
limits become less stringent as mA increases. 


• Antiproton flux bounds. Though B has baryon number +1, we nevertheless expect 
to obtain antiprotons from the AA —> Bcj) process, since the B decay products will 
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hadronize. 16 This additional antiproton flux can be tested in cosmic ray experiments 
like PAMELA [84] and AMS-02 [85]. The flux of antiprotons is given by [86] 

* (A> ® ) = S(S) (5.4) 

where K is the kinetic energy of the antiproton (a distinction important in the low 
energy limit), Vp is the velocity of the antiproton, and R(K ) is a best fit function that 
describes the propagation of the antiprotons throughout the galaxy (see [86]). As in 
the gamma ray case from Sec. 4, we can extract the antiproton spectrum dNp/dK from 
B decays in Pythia and boost to the AA rest frame. 

Various groups have derived antiproton bounds using AMS-02 data [87-89], typically 
showing results for WIMPs that annihilate to bottom quarks. To a reasonable approx¬ 
imation, the antiproton yield in B —> ebb decays is comparable to that of a bottom 
quark, yielding around 0.3 antiprotons per decay. More accurately, a single B —>• ebb 
decay from A A -A Bcj) yields a factor of 2.5 fewer antiprotons than an equivalent en¬ 
ergy bb pair. Thus, instead of evaluating Eq. (5.4) directly, we can simply scale down 
the XX bb bounds by this factor. Taking the “Ein MED” bounds from [87], for 
niA = 60 GeV we estimate 

(aAA^B<t>v) eS < 2.3 x 10 _26 cm 3 /sec, (5.5) 


which is in some tension with the GC best fit region. That said, there is at least a 
factor of two or three astrophysical uncertainty in these bounds. In addition, while the 
CMB bounds are irreducible, in the sense that the same photons that explain the GC 
excess will inevitably correspond to power injected into the early universe, perturbing 
the CMB, the antiproton bounds are less directly tied to the gamma ray signal. For 
both of these reasons, we have opted not to show the antiproton flux bounds in our 
plots. 


• Direct detection bounds. Because A has couplings to 4> and 4> couples to SM states, there 
will necessarily be a contribution to A-nucleus scattering from t-channel (f) exchange. 
For specific models, such bounds might be relevant (and prospects for future direct 
detection experiments promising), though direct detection constraints can typically 
be avoided for two reasons. First, while </> needs to have large enough coupling to 
stay in thermal equilibrium with the SM, those required couplings are small from the 
perspective of A-nucleus scattering. Second, small increases in can lead to large 
decreases in the A-nucleus cross section, since f-channel scattering typically scales like 
l/m^ at small recoil energies. 

For the specific model studied in App. A, <p mixes with the SM Higgs after electroweak 
symmetry breaking with a mixing angle O^h- The spin-independent scattering cross 
section of A through </> is (see, e.g. [57, 90]) 


osi - 




<f)h 


4 „,2 


7 r rrr, v 


(f) EW 


(5.6) 


16 Because of its baryon number, every B decay will necessarily lead to at least one proton (directly or from 
heavier baryon decay). Such a proton excess, however, does not seem to be visible over cosmic ray proton 
backgrounds. 
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where An is the coupling of A to cj>, dew is the Higgs vacuum expectation value, / is a 
factor obtained from the different parton fractions which we take to be 0.35 [91], and 
m n is the nucleon mass. Taking = 250 MeV as a benchmark, this cross section 
scales like 


<7Si ~ 9 x 10 46 cm 2 


/An \ 2 /250 MeV \ 4 

VHOy V ri H ) 


e, 


(f)h 


1.1 x 10 


-7 


2 


(5.7) 


where the baseline value saturates the LUX bound [92]. While this baseline mixing 
angle is small, <j) still decays prior to BBN. Of course, larger mixing angles are allowed 
for larger m^, though one should then account for (j> decays in the GC excess analysis 
(see footnote 14). 


• Flavor bounds. One aspect of this scenario that we have not delved into deeply is the 
generation of the B decay operators in Eq. (2.4), and in particular their flavor struc¬ 
ture. At the scale A (or below), some new heavy states are required to generate these 
operators, and those states could contribute to flavor-violating interactions. For any 
Bu c d c d c -like operator and assuming that A-scale physics conserves CP, one expects that 
the strongest limits should come from meson mixing induced by those new heavy states 
[93-95]. Since all fields involved in the B decay are right-handed, stronger constraints 
from chirality-mixing operators can be avoided. Assuming a generic flavor structure for 
the new physics, the most constraining bound comes from the (cr^ur ) 2 / A 2 operator, 
which is bounded by A > 1200 TeV [94], Taking ms = 60 GeV in Eq. (2.5), this 
yield a B decay width of T# ~ 1.2 x 10 ~ 20 GeV. Interestingly, T r ~ H occurs at a 
temperature of 121 MeV, above the beginning of BBN at T ~ 10 MeV [46, 61], so this 
flavor-safe scenario is indeed consistent cosmologically. Smaller B lifetimes (such as 
for the A = 300 TeV benchmark we use in Eq. (2.5)) require some mild suppression of 
flavor violation among the heavy states, which is certainly plausible if A-scale physics 
is approximately flavor conserving. 


• Collider searches for displaced jets. The operators that lead to B decay can also lead 
to B production at the colliders. For our B —> ebb benchmark, this cross section is 
rather small at the LHC, but if the B —> udd operator is present, then there is a more 
promising process: 

ud -» Bd. ( 6 - 8 ) 


For the flavor-safe case of A ~ 1200 TeV, this cross section is negligible, but being 
optimistic about flavor bounds (and pushing beyond the recommended values in Fig. 5), 
we take A ~ 30 TeV as a benchmark to explore possible LHC signatures. Since ms — 
60 GeV and B decays to quarks, this yields a relatively low energy four jet final state, 
and it is questionable if such events could be seen over overwhelming QCD backgrounds. 
That said, plugging A = 30 TeV into Eq. (2.5), the B has a lifetime of tr — 2.2 x 
10~ n sec, or a decay distance of ctr — 0.7 cm. Thus, the jets from B decays come 
from a (potentially very) displaced vertex, similar to the phenomenology of hidden 
valleys [96, 97] and their variants (see e.g. [98-100]). 

To get an estimate of the B production rate, the parton-level cross section in Eq. (5.8) 
scales like 


(7 (ud —> Bd) ~ 


1 


is — m 


2 \2 
B) 


167T 4sA 4 


~ 47 x 10 -2 fb 


(14 TeV ) 2 


30 TeV 
A 


(5.9) 
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To obtain the proton-proton cross section for the 14 TeV LHC, we integrate over the 
MSTW2008 LO parton distribution functions [101]: 



(5.10) 


For the high luminosity LHC, with target luminosity of 3 ab \ we expect around 120 


events. With a dedicated displaced jet trigger, one could hope to identify these events. 
This would be a distinctive signature for this scenario, especially if one could somehow 
verify that every B decay yields a baryon +1 final state. We leave a study of the LHC 
detection prospects to future work (see related studies in [102-104]). 

6 Conclusions 

In this paper, we presented a DM scenario where the interactions responsible for asymme¬ 
try sharing in the early universe are potentially visible today through indirect detection 
experiments. The key novelty compared to other asymmetric DM scenarios is that the DM 
abundance is set by thermal freeze-out involving an unstable particle B rather than by the 
decoupling of high-scale interactions. Assuming that B chemically decouples above the DM 
mass, the parametrics of A freeze-out behaves much like a standard WIMP, albeit with a 
non-zero chemical potential and a rescaled freeze-out value xp —> yp. Given this connection 
to WIMP physics, it is not surprising that we found a benchmark scenario with the right DM 
and baryon abundances. Intriguingly, this same benchmark yields a gamma spectrum from 
A A —>• Bcj) semi-annihilation compatible with the GC gamma ray excess seen by Fermi. 

This work emphasizes that non-minimal asymmetric DM scenarios can produce inter¬ 
esting indirect detection signals without relying on a residual symmetric DM component. 
As shown already in Refs. [41, 43, 44], asymmetric DM models with multiple stable states 
can yield indirect detection signals from semi-annihilation. Here, we emphasize the role that 
unstable dark sector states can play both in generating indirect detection signals as well as 
in sharing primordial asymmetries. 

An interesting variant to our scenario is if B were stable, perhaps due to an additional 
Z 2 symmetry. In that case, both A and B would contribute to the DM abundance, with the 
relative ratio determined by the A A —> B<p process. As long as there is a sufficient abundance 
of A particles, then this semi-annihilation process could give rise to both a boosted DM signal 
from the final state B [105-108] as well as a standard indirect detection signal from the decays 
of <fi. We leave a study of this scenario to future work. 

Finally, the nature of DM cannot be determined through any single observation, and 
even if the GC excess is indeed due to DM (semi-)annihilation, one would want to determine 
the dark sector properties through a suite of other experiments. For this particular DM 
scenario, the hadronic B decays imply an irreducible cosmic ray signal from antiprotons and 
potentially antideuterons. More intriguingly, evidence for B particles might show up at the 
LHC or future colliders through displaced hadronic decays, yielding a visible portal to a rich 
dark sector. 
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Field 

Spin 

Baryon number 

ll>A 

Weyl left 

+ 1/2 

r A 

Weyl left 

- 1/2 

4>a 

complex scalar 

+ 1/2 

<t>°A 

complex scalar 

- 1/2 

B 

Weyl left 

+1 

B c 

Weyl left 

-1 


complex scalar 

0 


Table 1. Field content for the concrete model, which generate the required interactions specified in 
Figs. 3 and 4. 
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A Concrete Model 


For simplicity in the text, we referred to A as being a single particle. Since the decay 
operators in Eq. (2.4) imply that B must be fermion, though, A should be replaced by a 
fermion/boson system of iI’a/ 1 j>A- Motivated partly by supersymmetry but mainly by the 
resulting simplification of the Boltzmann equations, we take Vm/V’a to be a Dirac fermion 
that is mass degenerate along with two complex scalars 4 > a/4> c a (i.e. two chiral multiplets 
with a holomorphic mass). We comment on the impact of splitting the Vt 4 /(ftA masses below. 

In the text, we also referred to 4> as being a single particle and used the benchmark 
m$ = 250 MeV for studying CMB and direct detection bounds. Here we replace 4> with two 
real fields, one scalar and one pseudoscalar. This can also be motivated by supersymmetry, 
but more relevant to our scenario, if cj) were composed of a single real scalar field, then the 
annihilation i/ja^A —3► would be p -wave suppressed. We consider these two fields to have 
somewhat heavier masses than considered in the text, in order to avoid additional bound 
state formation (see footnote 7). We will see that this has a negligible effect on the CMB 
bounds in Eq. (5.3) and weakens slightly the direct detection bounds in Eq. (5.6). 

The field content for this model is shown in Table 1 , where B/B c is a Dirac fermion and 
cj) is a complex scalar. We separate <p into its scalar component s and pseudoscalar component 

(2, 8/S 


s + ia 
s = — 


(A.l) 


This is similar to the axion portal [55, 56], though we have suppressed a possible vacuum 
expectation value (vev) for <p in order to remain agnostic as to whether or not a is a pseudo- 
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Goldstone boson from spontaneous symmetry breaking. The Lagrangian is then given by 


free 


r 4 

^mt 


r 6 

^mt 


iip^d^A + i^ C A a^d fl ^ c A - (m A ip c A ^A + h.c.) 

+ |^ A | 2 + |^| 2 - mi (|^| 2 + |^| 2 ) 

+ iBo^d^B + iWa^d^B 0 - ( m B BB c + h.c.) 

i a m2 1 2 2 1 2 2 

+ Wl “ 2 m ' S ~ 2 m “ a ’ 

A ab(iPaB c (Pa + ^aB^a + h.c.) 

- (A A'lpA'ipA^ + A bBB c ()) + h.c.) 

+ ^A<j> (|0a| 2 + \<Pa\ 2 ) (4> + $) + A A<t> (|^a| 2 + I^aI 2 ) $<P 
Bu c d c d c , 


(A.2) 


(A.3) 

(A.4) 


Note that we have introduced a mass splitting between the scalar and pseudoscalar modes, 
such that the decay s —> aa is kinematically allowed. For simplicity, we have taken the 
Yukawa couplings to s and a to be the same, and assumed various relations among the scalar 
4 >A and 4 >a couplings. 

The potential term V ( H , <p) includes terms that mix cp with the Higgs sector after 
electroweak symmetry breaking. Motivated by supersymmetry and the axion portal, we 
assume a two Higgs doublet model with H u d , such that one possible mixing term is 


V ( H , q i>) D fJtpH 4> H u H d + h.c. 


(A.5) 


This interaction assures that both s and a can decay to SM states, though V(H, cf) generically 
contains cubic interactions such that the s —> aa decay mode dominates. We ignore the 
induced vev of 4 > from this mixing, since it can be absorbed into a redefinition of the other 
couplings. Depending on the potential, the Higgs vevs can also contribute to the </> mass, 
though we expect this to be a small effect given the small mixing angle suggested by Eq. (5.6). 
Without supersymmetry, some degree of fine tuning would be necessary to keep the s and a 
masses small given their large couplings to the A and B states. 

The baryon assignment of the different particles is set by Bu c d c d c / A 2 in Eq. (A.4). So 
unlike in the SM, baryon number is not an accidental symmetry of this Lagrangian, though 
that conclusion might change depending on the precise dynamics present at A. 

For the purposes of App. D, the key feature of the coupling choices above is that the 
A scalars are treated symmetrically, such that any process involving c/)a has a counterpart 
involving (^ 4 ) ■ This, along with the assumed mass degeneracy of the A fermions and A 
scalars, allows the full Boltzmann system to simplify to a two-particle system. We checked 
that all the relevant DM interactions from this Lagrangian include an s-wave term (see 
e.g. [109, 110]), which is necessary to achieve the desired cosmology. For example, an s-wave 
annihilation channel for the A fermions is possible via 4’a' 1 Pa sa. 

Introducing mass (or coupling) splittings between ip a an d <pA would add new parameters 
to the model studied and complicate the Boltzmann equations. The dominant effect of such 
a splitting is to change the resulting populations of the ip a and <p A components. To maintain 
roughly the same phenomenology as presented in the text, we need m,/ M /rn ( j )A to be 0 ( 1 ), 
such that the lighter DM component is still Boltzmann suppressed at freezeout. We also 
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need to satisfy the mass ordering 

max{m^,m^} < m B < m^ A + m ( j >A . (A.6) 

This first inequality comes from the requirement of not depleting A particles by scattering 
with B particles prior to freezeout (see footnote 13), and the second inequality ensures that 
the semi-annihilation process i^A^A H► Bcj) is kinematically allowed even at threshold. 

In addition to the possibility that B might correspond to a bound state of A particles, 
the 4> mediator might give rise to additional bound states of A and/or B. For simplicity of our 
analysis, we want to avoid the possibility of />-mediated bound states, which would further 
complicate the Boltzmann analysis. Note that only the scalar s mediates a 1/r Yukawa 
potential between the A particles, whereas the pseudoscalar a can only mediate a 1/r 3 spin- 
dependent potential among the A fermions. For our benchmark, with 0(1) (j) couplings and 
DM masses around 50 GeV, we can avoid bound state formation if m s > 0(1 GeV). 

The presence of both s and a turns out to have relatively little effect on the phenomenol¬ 
ogy presented in the body of the paper. There are now two semi-annihilation processes rel¬ 
evant for indirect detection, i/ja^a —> Ba and il’A^A -A- Bs. As long as 2< m a < 3m n 
(such that a dominantly decays asa-> H + H~) and m s > 2 m a (such that s dominantly decays 
as s —> aa —> h + H~ H + H~), then the resulting photons from muon FSR give subdominant 
contributions to the photon flux as needed for Sec. 4. The CMB bounds in Eq. (5.3) only 
have a mild dependence on the injected muon spectrum so are largely unchanged. Raising 
m s weakens the spin-independent direct detection bound in Eq. (5.6) and spin-dependent 
bounds from a exchange are typically sub dominant. 

Finally, we note that similar phenomenology could be achieved by replacing with a 
U(l) / gauge field that kinetically mixes with SM hypercharge [48-50]. The main challenge for 
using this hypercharge portal is that the B decay operator in Eq. (A.4) would then require 
an insertion of the U(1/-breaking Higgs field, making it an even higher dimension operator. 

B Chemical Potential Analysis 

To determine the abundance of DM at the decoupling temperature To, we have to consider 
the chemical potentials of all relevant SM and dark species. For simplicity, we work in a 
regime where the sphaleron process becomes inactive prior to To, such that we can treat the 
baryon and lepton asymmetries as being independently conserved during B decoupling. The 
states in the dark sector then carry effective baryon number consistent with the operators in 
Eq. (2.4). 

In the early universe, SM interactions guarantee chemical equilibrium among SM par¬ 
ticles. After electroweak symmetry breaking, the chemical potential relations can be written 
as [62, 111] 

Hu Hd = Hu — He = HW i (H.l) 

3 (Hu T Hd) — Hu 4“ He, (B.2) 

where Hu, Hd, Hu, and He refer to the chemical potentials of each flavor of up-type quark (ul 
and ur), down-type quark (cLl and cLr), left-handed neutrino (vl), and charged lepton (e^ 
and cr). Note that Eq. (B.l) is enforced by exchanges (with chemical potential hw) 
and Eq. (B.2) is imposed by the sphaleron process. At temperatures below the top quark 
mass, the neutrality condition imposes 

8 Hu ~ ®Hd - 6 He + 6 Hw = 0. (B.3) 
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The above relations imply 


M u — ^Md — ig Mej (B-4) 

which allows us to write the chemical potential for SM baryons as 

90 

M& = (3 — l)(Mui, T Mur) + 3(Mc/l + Mc!r) = ~ JgMe j (B.5) 

where we do not include the contribution from the top quark. The semi-annihilation process 
in Eq. (1.1) and the transfer operators in Eq. (2.4) impose the relations 

MB = 2/m, (B.6 ) 

Mb — M« T 2 fid- (B.7) 


Using Eqs. (B.4) and (B.5), the chemical potentials for A and B can be simplified to 

29 

MB = 2/m = — /M 


(B.8) 


For a species i with gi degrees of freedom and mass m;, the relation between the charge 
density and the chemical potential /q at temperature T is [112] (see also [40]) 


ni - n ■ = g t f(rm/T ) T 3 


(B.9) 


where we are assuming /q -C T with 

1 TOO 

47t 2 Jm,' 


V 2 dy 


\ _ J 4?r2 Jmi / T cosh 2 (\^/x 2 +y 2 ) 

J\ J ') ^ 1 roo _ y 2 dy _ 

4iU Jm,i/T g^h 2 


= /f(x) (for fermions), 
= /b(x') (for bosons). 


(B.10) 


The function f(x) takes into account the Boltzmann suppression of particle i. In the limit 
x 2> 1, f{x) ~ 2(x/27r) 3 / 2 e _a: . Using Eqs. (B.5)-(B.10), we can write the asymmetries as 
r}i = Xi(m — rij)/s, where Xj denotes the baryon charge of species *. This leads to 


Vb 

MM 

W 

MB 


1 T 2 /i 6 


/f(0) 


1 29 g^ A T 2 g b 


2180 s 
1 29 9 ^ T 2 g, b 
2180 1 
29 gs T 2 g b 
90 s 


fi{m A /T ), 


- Mm A /T) 


ff(m B /T ), 


(B. 11) 
(B.12) 

(B.13) 

(B.14) 


where g q is the number of degrees of freedom for each quark flavor, i.e. g q = 3 from color (since 
the left- and right-helicities are already included in Eq. (B.5)). For the system described in 
App. A, 

MM = 9</>a = 9<f>\ = 1) 9b = 2. (B.15) 

The relations above are valid for temperatures T > T B , where is the B decoupling 
temperature. Using 

Va = MM + Mm + W A , (B. 16) 
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the total asymmetry r/tot in Eq. (3.1) can be written as 


f?tot = VA + VB + Vb 

[ 29 3 (g^ A ffjmA/T) + 2g (j)A f h (m A /T)) 29 3g B f{(m B /T) ' 
~ T]b UeO g q f f (0) 90 g q ft(0) 

(B.17) 


which is the basis for Eqs. (3.2)-(3.4). From Eq. (B.9) and the relation above, it is also 
possible to write the DM chemical potentials as functions of rjtot as 


Mi _ f 27t 2 g* s \ Vi 
T ~ V 45 gi ) f(rm/T ) ’ 


(B.18) 


where we used s = 2ir 2 g its T 3 /45. 


C Thermal Cross Section for Bq —> qq 


The process Bq -A qq , where q is a quark consistent with the operators in Eq. (2.4), is efficient 
at high temperatures, and must be included when determining the B decoupling temperature 
To- From the amplitude 

— 'i 

iM = -p (zix 2 )(y32/4) (C.l) 

expressed in terms of the Weyl wavefunctions x, and yi, we find 


1 

4 




s(s — m 2 B ) 
4A4 ’ 


(C.2) 


where s is the square of the center-of-mass energy, and we assume mg> m q . The differential 
cross section with respect to the Mandelstam t variable is 


1 


\Mf 


da 

dt 167T A(s, m 2 B , 0) ’ 

where to simplify the phase space factor we use the function 

A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2 xz — 2yz. 
Integrating the cross section over t, we find 

s 

a = 


647tA 4 


(C.3) 


(C.4) 


(C.5) 


At sufficiently high temperatures, we can ignore the chemical potentials of B and the 
quarks. For a generic scattering process 12 —> 34, the Mpller thermal cross section is [113] 


(crDviol) 


1 

12—>-34 cq oq 

n 1 n 2 


9l dV g 2 d 3 p 2 (E 1+ E 2 )/T 

(2vr)3 (2 tt) 3 Mo1 


(C.6) 


where the Mpller velocity is defined such that the product umo! 714712 is Lorentz invariant, and 


n?(T) = 


_ f 9i d 3 Pi -e/t 


(2vr) 


(C.7) 
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(C.8) 


In our case, n^ q = n e B0 is the thermal Boltzmann distribution for B , 


eq _ 
bo ~ 


^(^VW^T 3 , 


2tt 2 V T ) 


\ T > 


where Ki is the modified Bessel function of order i, and re 2 = is the relativistic distribution 
for quarks, 

n q0 = ^ C ^g q T \ (C.9) 

where here we have corrected for the Fermi statistics of the relativistic quarks. The rela¬ 
tive velocity UMoi is obtained from p\ • P 2 = 'UMoi-Ei-E^. The integral in Eq. (C.6) can be 
reparameterized following the annihilation example of Ref. [113] as 


d 3 pi d 3 p 2 = 2tt 2 E\E2 dE + d E- ds, 
with the limits of integration: 

Ejf. = E\ + E 2 > y/s, 

E- = E\ — E 2 set by the condition | cos0| < 1, 
s > m 2 B . 

Performing these integrals, (ctumoi) Bq^qq is given by 

<™m„i )B-^ m = ^ ^ 32 A 4 f m% is ^ s - m *? K 1 (t) 

48vrA 4 C(3)^ 2 (®a) ' 


(C.10) 


(C.ll) 

(C.12) 

(C.13) 


(C.14) 


D Boltzmann Equation Details 

In going from the full model in App. A to the simplified discussion in Sec. 2, we asserted that 
the four-particle system X £ {il>A, 4 >a, (t>% B} could be reduced to an effective two-particle 
system. Indeed, this is possible if all states type A are mass degenerate and have equal 
chemical potentials such that p^ A = = ha- We now derive the full Boltzmann 

system and show how this reduction occurs. 

Unlike the discussion in Sec. 3.3 where we introduced the y variable, here we stick to 
the notation x = ra^/T. The equilibrium distributions for a state X are given by 

Y x^ x ) = Y xo( x ) ex P(/u7 T ), ( D -l) 

where yx is the chemical potential of X and we assume yx *C T which is a valid assumption 
until freezeout. In the non-relativistic limit (mx T), bosons and fermions follow the same 
distribution [47], and the equilibrium functions denoted with a subscript “0” are 

= (D.2) 

We use the notation A = s/H(m,A), where s is the entropy of the universe, H the Hubble 
constant at the temperature T = tua, and g* s is the effective number of relativistic degrees 
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of freedom at the same temperature [47]. The chemical potentials satisfy 2 //a = Hb when 
the semi-annihilation process AA — > Bcj) is in equilibrium. We also have //a = —/i 4 and 
Hb = —/Cg when the annihilation processes AA — > 00 and BB —> 00 are in equilibrium. 
Here we assume = 0, since the 0 is in thermal equilibrium with the SM and there is no 
symmetry forbidding interactions like 00 —> 000. 

By assumption in App. A, the couplings involving 0a and (jfj^ are identical, so we know 

that 

U, = y k (d. 

at all temperatures T. Therefore, processes that involve 0a, such as s i m ply 

get a factor of 2 from including the equivalent process with 0]J. With this simplification, the 
Boltzmann equations for the {0 a, B} system are: 


dx 


d Y b 

dx 


x* 




( Y e q y e q \ 

Y*aY* a -Yb^£± j 


/ y e q y e q x 

- 2 W^b^a^ ~ Y 'i’ A a J 

/ y eq yeq \ 

a .V) l t -%#- I 


+ 2 (&. 




i’A 1 B 1 6\ y e q J 

1 / J 




6 


yeq yeq > 

- 2 {vm^b+v) I Y^ a - Yb^P^ 


( yeq yeq >. 


yeq yeq 

A B 


\ 


+ 2 A^ a B^ a A S> \ V A Y B y ^ a y *[ J 


+ ( cr Bg^5g v )^SM ( Y B ¥ B q ) 


T B X 

H 


(Yb ~ Y^) 


(D.4) 


(D.5) 


The Boltzmann equations for 0 a can be easily extrapolated from those of 0 a- Here, we have 
neglected the conversion process A A —> BB since annihilation dominates over conversions 
for determining the freezeout dynamics of near-mass states (recall that our mass range of 
interest is tha < ms < ^tua)- 

At high enough temperatures (though not so high that boson/fermion statistics matter), 
the equilibrium distributions satisfy 


V eq _ 


2Y Z 


(D.6) 


where the factor of 2 comes from the number of degrees of freedom. In general, though, the 
evolution of Y^ A and 2 Y0 4 will not remain identical, except for special choices of the couplings 
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in Eq. (A.3). To see which combination is needed, consider the evolution of Y^ A —Y ( p A — Y^t = 

Y^ a ~ 2 Y^ a . 


dx ^ 2Y * a) ~ x 2 




i’A 


<t>A<t>A^ r <t > 4 1 


+ 2 (% 5 *») («* - 2 ^) 1 % - (yI - 2 xly ?< ) 


— 2 (a A r> . , .u) 




y :yt\ 

yeq I 

<t> ] A ' 

yeq yeq 
B 


v e q 

^A 


(D.7) 


In addition to the factors of 2 coming from including (p A -H- </)j, there is an additional factor 
of 2 in the second line coming from averaging over polarizations in the cross section, since 


W A B -+^> = 2 (% A B^U V) - (D ‘ 8) 

There is a factor of 1/4 in the second line from the ratios of the equilibrium distributions 
since = 21/l/\ In the third line, the factor of 1/2 comes from averaging over spins in the 
cross sections, and the factor of 4 comes from the ratio of equilibrium densities. From the 
semi-annihilation terms alone (i.e. second and third lines), it is now clear that if Y ^j = 2 Y^ 
at early times, then this relation will continue to hold at later times. For the annihilation 
terms (i.e. first line) to maintain the same relation, we simply need to choose couplings such 
that 

= 2 KaW 2^4> v ) ■ ( D - 9 ) 

This in turn sets the needed couplings in Eq. (A.3). Of course, for more general couplings, 
one can simply solve the full set of Boltzmann equations. 

Assuming Eq. (D.9) holds and plugging the relation Y^ A = 2 Ya a into the Boltzmann 
equations, we arrive at an effective two-particle system with {A, B}: 


dY A 

dx 


d Yb 
dx 




(iT ) 2 


T {&ipA<t>A —( Y A ^ B ycq 


- Y a Y b ~ Y A 




B 

Y?Y?\ 

YT J 

ye q y eq ^ 


B ± A yeq 




(oBB^v){Y B Y-s-Y?Yg) 

Wif>A<l > A—^Bif> v ) — Yb yeq ^ 


(D.10) 
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y_eqy eq \ 

A B \ 


+ 2 P y * Yb - Y * wy I 

+ {a m ^,v)YZ (Yb - Y?) - T -f- (Yb - Yj>) . 


(D-ll) 


with conjugate equations for A and B. Here we refer to Y^ A as simply Ya for simplicity, but 
we maintain old notation for the cross section in order to maintain the correct factors of 2 
in averaging over polarizations. 

As one can check, the B decay term and the Bq —> qq scattering term depend explicitly 
on the chemical potentials, as expected from the discussion in Sec. 3.2. Below Tp>, though, 
when the scattering term can be neglected, the chemical potential factors cancel in the 
Boltzmann system, since we take ratios of the equilibrium densities except for the decay 
term. In other words, after Tjj, we could make the replacement Yj q —>• Y ^-q, as defined in 
Eq. (D.2). Numerically this is equivalent to setting the chemical potentials solely as initial 
conditions, as given by, say, Eq. (3.8). 
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